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CO . Abstract 

Physics interpretation of the "magic" baseline, L mag i c , that can play important 
qq | role in future oscillation experiments is given. The "magic" baseline coincides with 

the refraction length, Iq. The latter, in turn, approximately equals the oscillation 
length in matter at high energies. Therefore at the baseline L = Iq the oscillation 
phase is 2ir, and consequently, the "solar" amplitude of oscillations driven by #12 and 
Am|i vanishes. As a result, in the lowest order (i) the interference of amplitudes in 
the v e — v„ (v T ) transition probability is absent; (ii) dependence of the probability on 
the CP-phase, 5, as well as on 9 12 and Amfi disappears. Corrections to the equality 
1 1 Lmagic = lo are estimated. Effect of changing density is considered and two new magic 

^h. trajectories are identified for neutrinos that cross the core of the Earth. Other magic 

baselines associated with zeros of the atmospheric amplitude are discussed. 
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1 Introduction 



It was observed some time ago that at the baseline 

ma9lc ~ V2G F n e } [ ) 

where Gp is the Fermi coupling constant and n e is the electron number density, the analytic 
formula for the — v e oscillation probability in matter (in 3z/— mixing context) takes a very 
simple form J]Q. The probability does not depend on the CP violation phase, 5, as well as on 
the mixing angle, 9u, and the mass splitting, Am| 1; of the 1-2 sector. Therefore, neutrino 
oscillation experiments with the baseline L = L mag i C = (7300 — 7600) km will allow one to 
perform clean measurements of #13 resolving degeneracy with the phase S PUI2]- That was 
discussed in a number of recent recent publications in some details The baseline L magic 
depends only on matter density and does not depend on neutrino energy and oscillation 
parameters. For this reason it was termed the "magic" baseline in 

Apparently the "magic" baseline defined in (0) coincides with the refraction length, Iq, 
introduced by Wolfenstein almost 30 years ago [3]: 

Lmagic = ^0- (2) 

This is not accidental coincidence and in what follows we will explain the reason behind the 
equality (0). The explanation is given in sees. 2 and 3. Sees. 4 and 5 contain some more 
advanced material, in particular, discussion of various corrections to the equality (J2J). 



2 Magic baseline and refraction length 

Recall that the refraction length, l , has been defined as the distance over which an additional 
phase difference, (f> m , acquired by neutrinos due to interactions with matter, equals 2n: 

m = VIq = 2tt. (3) 

Here V = \f2Gpn e is the difference of potentials for v e and in usual matter, so that Iq is 
a characteristic of medium relevant for v e — mixing. 

The refraction length enters the MSW- resonance condition [B], 

Z o cos20 = Z„, (4) 

where l v = 4nE/Am 2 is the vacuum oscillation length and 9 is the vacuum mixing angle. 
Inverse quantity, 1/Zo, determines the eigenfrequency of medium, and for small mixing the 
condition (jjj) means that the eigenfrequency of medium coincides with the eigenfrequency 
of neutrino system, \jl v . For large mixing (strongly coupled system) there is a shift of 
resonance frequency given by cos 29. 
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The ratio Iq/I v determines modifications of the mixing angle, 9 m , and oscillation length 
in matter l m : 

sin2# m = sin2#[(cos2fl - l v /l Q f + sin 2 2fl]- 1/2 . (5) 

l m = q(cos20 - l u /l ) 2 + sin 2 26}- 1/2 . (6) 

According to (jfi|). at low energies (below resonance) one has l m ~ l v \ with increase of energy 
the length l m first increases, it reaches maximal value, l™ ax = lo/ sin 29, at l v = Iq/ cos 29 
(i.e., above the resonance: E max = Ep/ cos 2 29) 1 , and then decreases approaching Z from 
above. In the non-resonance channel, l m increases with E approaching lo from below. Thus, 
in the case of large densities of matter or large energies of neutrinos, when l v ^> l , we 
obtain 

lm « lo- (7) 

In other words, in the matter dominating case the oscillation length in matter approximately 
equals the refraction length. This is the key point of physics interpretation of the magic 
baseline in the next section. 



3 Where does the magic baseline come from? 

Let us consider oscillations of three mixed neutrinos Vj = (i/ e , z/ M , z/ T ) T in the matter. The 
vacuum mixing matrix, Upmns-, that relates v$ and the mass eigenstates v = (ui, u 2 , vz) T 
Vf = Upmnsv, can be parametrized as 

UpMNS = U23hUi3l-sUi2- (8) 

Here Uij = Z7y (%) performs rotation in the ij- plane by the angle Oij and Is = diag(l, 1, e lS ) 
is the matrix of CP-violation phase. 

The —>■ v e transition probability can be represented as 

P{v» -^v e ) = \ cos 9 23 A s e lS + sin 9 23 A A \ 2 (9) 

(see some details in sec. 4). Here As is the solar amplitude that in the lowest order 
approximation (up to corrections of the order Am^/Am^, and sin 2 #i 3 ) depends on the 
solar neutrino oscillation parameters: Am^, #12- A a is the atmospheric amplitude that 
depends on the atmospheric neutrino oscillation parameters, Am^, #13. The CP- violation 
effects are given by interference of the two amplitudes in Q. 

Let us consider the constant density medium that, in fact, is a very good approximation 
for neutrinos propagating in the mantle of the Earth. Up to the phase factor we obtain for 
the uniform medium 

A 5 = sin 20™ sin (10) 

1 I am grateful to J. Kersten for correcting the maximal value. 
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where O 7 ^ is the 1-2 mixing angle (0) and (j) m is the oscillation phase in matter: 

M - ^; (ii) 

I'm 

the oscillation length, Z m , is given in (JBJ). Apparently As is a square root of the usual 
oscillation probability. (Similar expression can be written for Aa with substitution 1-2 
parameters by 1-3 oscillation parameters.) 

Let us consider neutrinos with energies E > (0.5 — 1) GeV relevant for accelerator exper- 
iments and atmospheric neutrino studies. For these energies ^(Am^) = AnE/Am^ ^> lo, 
and therefore the solar oscillation mode is in the matter dominating regime when 

l m {Am 2 21 ) « lo- (12) 
Correspondingly, the phase of oscillations equals 

W - (13) 







and at the baseline L = l we obtain 



2vr, A s = 0. (14) 
The solar amplitude vanishes, and consequently, 

P{y^v e ) « \sm8 23 A A \ 2 . (15) 

So, for high energies and L = l the "magic" properties are reproduced: 

- the interference and dependence on the phase 6 disappear; 

- the dependence on solar oscillation parameters disappear too. 

Summarizing, the "magic" length is nothing but the refraction length, at least in the 
lowest order in small parameters. For high energies at the distance L = l Q , the phase of "so- 
lar" oscillation amplitude becomes 2n. The amplitude vanishes, and therefore dependence 
of probability on CP-phase disappears. Similar consideration is valid for the v e — v T channel 
with substitution sin 6*23 — > cos #23 and cos 6*23 — > — sin #23 in Q. 



4 Corrections to the equality L mag i c = Iq 

In the previous section we have neglected terms of the order sin 2 613 and Am^/Am^ which 
can be as large as 3%. Also we have taken l m = l , but at relatively low energies the devi- 
ation from this equality can be significant. Notice that the 3% correction to the refraction 
length of 7300 km equals 220 km which may be non-negligible for selection of a detector 
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place, if high precision measurements of parameters are planned. One way to proceed is to 
perform numerical search of the baseline when the probability has the weakest dependence 
on S j2j. Here we present some analytic consideration. 



We will define the magic baseline as the distance over which the "solar" amplitude 
vanishes: 

As(L mag i C ) = 0. (16) 

In the constant density approximation (and in the lowest order in sin #13) this leads to the 
condition L magic = l m , and for high energies l m ~ l . There are two types of corrections to 
the equality L magic rs l : (i) due deviation of l m from the asymptotic value, l , and (ii) due 
to dependence of As on sin# 13 . We consider these corrections in order. 

1). Corrections due to l m 7^ 1$: For l v l we have 

Lmagic = kn « l (l + COS2012~) = l (l + COS29 ^^^) ( 17 ) 

in the resonance channel. The correction is positive and equals ~ 0.1/ E (GeV). For a con- 
stant density medium, with increase of energy the correction and the magic baseline become 
smaller. However, in the case of realistic Earth density profile an average density along a 
trajectory increases with a length of trajectory, or equivalently, with | cos where 9^ 
is the zenith angle of the trajectory. Therefore Iq decreases. So, for the Earth the magic 
trajectory does not change significantly with the energy being at | cos0j,| ~ 0.6. For the 
non- resonance channel the correction is negative (minus sign in Eq. (J17j0 . 



2). Corrections due to 1-3 mixing. To evaluate these corrections we need to give precise 
definitions of the amplitudes A s and Aa in Q and find their dependence on the oscillation 
parameters. The evolution of the flavor states Vf is described by the Hamiltonian 

H = UpMN S Diag(0, A 21 , A 31 )U j PMNS + V, (18) 

where V = diag(V, 0, 0), A21 = (S.m^^j'lE and A31 = Am| 1 /2£'. Let us define the neutrino 
propagation basis, v = (u e , z/ 2 , ^3) T , through 

Uf = U 2 zhv. (19) 

The Hamiltonian H that describes oscillations of v can be obtained from ([180. (JHJ) and (fT^j) : 

H = U 13 U 12 Diag(0,A 21 ,A 31 )Ul 2 U} 3 + V, 

or explicitly, 

( C 13 S 12^12 + V + sj 3 A 13 C13S12C12A12 S13C13A13 — S I3C13S-L2 A12 \ 
c? 2 A 12 -Si3Si 2 c 12 A 12 , (20) 

c 13^13 + S 13 S 12^12 / 
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where C13 
basis as 



cos 6*i3, S13 = sin ^13, etc.. Introducing the evolution matrix in the propagation 



S 



I A ■ -I I 



i,j = e,2,3, 



(21) 



it is straightforward to find that the — v e transition probability has the form with 



Aa = A 



e2- 



A A = A 



e3- 



(22) 



That is, the solar amplitude is given by the amplitude of transition i) 2 — > v e and the 
atmospheric amplitude coincides with the z/3 — > u e transition amplitude. In the limit 6*13 — > 
the state z/3 becomes the mass eigenstate z/3 and decouples from the rest of neutrino system, 
whereas z>2 becomes the combination of and v T that mixes with z/ e with the solar oscillation 
parameters. Therefore in this limit A e2 = A e2 (Am 21 , #12) as we discussed in sec. El 

For energies below the 1-3 resonance (E ~ 1 GeV), when the 33-element of the Hamil- 
tonian dominates, the corrections can be calculated immediately. Performing an additional 
rotation of the neutrino basis by ^13(^13) and then making block-diagonalization of the ob- 
tained Hamiltonian, we obtain that the heaviest state decouples and the two others form a 
2v— system with the mixing angle #12 and modified potential: 



V -> V 



s 2 13 V 



V 



A 
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Correspondingly, the magic length is modified as 



T ■ ~ 

-'-'magic ~ 



l + cos2# 12 ^ + s 



(23) 



(24) 



Notice that the block-diagonalization removes the imaginary part of the corrections to the 
solar amplitude (see discussion below). 



Let us present an estimation of the si3-corrections that are valid in whole the energy 
range including the 1-3 resonance, and also take into account the imaginary part of the solar 
amplitude. We perform an additional 1-3 rotation of the neutrino basis that vanishes the 
1-3 element of the Hamiltonian (|2(J|) : 

v = U 13 (e? 3 )u m . (25) 
Here v m = (^i m , i>2, ^3m)> an d the angle is given by 

tan 2^3= - 2He3 ~ , (26) 

with Hij being the r/'-element of the Hamiltonian ()20|) . In the new basis, the Hamiltonian 
has the form 

H£ cos(03 - 6 13 )a 12 
H m = I ... H 22 sm(6^-e 13 )a 12 j . (27) 
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where 



a 12 = s 12Cl2^ 



12- 



:28i 



Tjm 
#11 



cos 2 0%H„. - sm29[ n 3 H e3 + sin 2 6™H 33 , 



,2 nm 



H£ = sin^ 9? 3 H ee + sin 29? 3 H e3 + cos 2 9™H 33 . 



(29) 



Since we are looking for the corrections due to the 1-3 mixing the terms oc A12 can be 
omitted in the diagonal elements, H™, , and consequently we obtain 

. cos #i3 sin 0Jg 



Tjm 



COS #13 COS 9^ 



COS 



9 m 
13 



9 



13 J 



Tjm 



V- 



sm 



l™-9 



13, 



(30) 



Here we used expression for A13 in terms of mixing angle in matter that can be obtained 
from (EE 



A 



13 



1/ sin 20™ 



sin 2(0™ 
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13 J 



(31) 



In the propagation basis, the transition v e —>■ v 2 proceeds in two different ways: u e — ► 
v \m v 2 and u e —>■ v 3m — > v 2 . Therefore the "solar" amplitude can be written as 

A e2 = cos + sin (32) 

where A™ 2 and are the amplitudes of v\ m — > v 2 and z/ 3m — > v 2 transitions correspondingly. 
Then in the lowest approximation using the Hamiltonian (J27j) we obtain for (|32)1 : 



.4 



e2 



2a\ 2 cos ^ cos( 
+ sin 9™ 3 sin 



'13 







13 



-"11 



'13 







13; 



,sin(tf£L/2) _ <flJSL 



/2 



-"33 



(33) 



For energies much below the 1-3 resonance we have sinf?™, ~ sin #13 <C 1. Therefore the 
second term in (J3~3*)) is strongly suppressed, furthermore H™ ~ if ee ~ V, and consequently, 
the amplitude is reduced to the one considered in sec. 3 is recovered. For energies above 
the 1-3 resonance, cos 9™ 3 —>■ and the first term in (J32)) vanishes. Now « H ee « V and 
again we recover the result of sec. 3. 

Minimal value of A e2 corresponds to H^L w flgjijL ~ 27r when both terms in (J33)) are 
close to zero. (We confirm this by explicit calculation.) Then introducing small quantities 

1 2 11 > 3 2 

(ti <C 7r), and taking the first terms of expansions of sines and exponents in (}3*3*j) around tt 
we obtain: 

L 



"^33 L 



7T 



(34) 



(2 



2a 



12 



cosf?™. cos( 



+ sin#™sin(# m .-# 13 )(l- 
6 



I6ij 



7T 



rrm 
-"11 



^3) 



(35) 



The imaginary part contains an additional power of small parameters £j. 

Using explicit expressions for the H™ and (|30|) we obtain from (|35|) the real and 
imaginary parts of the amplitude: 

A? = Ci3Si 2 Ci 2 A 12 - y^j > ( 36 ) 



= -^Ci 3 S 12 C 12 A 12 V 

Let us analyze these results. 

1). The real part of the amplitude vanishes if 

2tt 



2tt\ 2 4vr 2 



+ ^ tan 2 13 



V V 2 



(37) 



Vci 3 (38) 

that coincides with the baseline obtained in ()24j) . 

2) . The imaginary part is always non-zero with minimum 

I^L^Imm = 27r 2 S 2 3 S 12 Ci2-^ (39) 

at L = 2-k/V. 

3) . At the baseline that corresponds to zero real part (|38jl . we obtain from (J37|) that 
the correction to the minimal value (f3T)|) is of the order sf 3 : [A^l = | m i n (l + s 2 3 ). 
This means that the baseline ()38)1 provides a minimum of the total amplitude in the order 
s 2 3 , and therefore it can be identified with the magic length. When corrections due to 1-3 
mixing are included, the solar amplitude does not vanish exactly, and therefore the magic 
properties are satisfied only approximately. 

4) . The amplitudes ()36ll3T|) do not depend on the mixing angle in matter, and therefore 
the results are valid in whole energy range including the 1-3 resonance region and the region 
above the resonance. 



5 The case of non-constant density 

In the case of non-constant potential (density) along the neutrino trajectory, the refraction 
length can be defined by the condition 

lo 

dxV(x) = 2tt. (40) 
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In the lowest order, for the "solar" amplitude in the matter dominating case one can obtain 
the following expression (see [7] for details) 



A s = Jsin2# 12 A 12 / cfeexp ( -i I dyV{y) ) . (41) 



2 

Consequently, the magic baseline can be found from the condition 

dxexpi^ij dyV(y))=0. (42) 

The double integration takes into account the change of both oscillation length and mixing 
angle along the trajectory. 

Let us consider the three-layer profile with constant potentials (densities) V m , V c and 
V m , and baselines L m , L c and L m . The densities and baselines in the first and the third 
layers coincide. To a good approximation that corresponds to profile along the neutrino 
trajectories that cross the core of the earth with V m and V c being the potentials in the 
mantle and the core correspondingly. 

Performing integration in eq. (J41|) we find 



As = sin 29. 



in 



sin| f + sin y 



(43) 



where 

sin20™«sin20 12 — - (44) 
is the mixing angle in the mantle, and 

4> m = V m L m , <f) c = V C L C (45) 

are the phases acquired in the mantle and the core. The sum C + 20 m is the total oscillation 
phase. The first term in f|^H|) corresponds to the amplitude of pure adiabatic transition. It 
depends on the mixing angle at the surface of the earth and on the total phase. The second 
term is the correction due to the adiabaticity violation at the border between the mantle and 
the core. So, in the presence of the adiabaticity violation the amplitude is not determined 
by total phase. 

According to ()43)1 the magic baseline is determined by the condition 

For constant density along whole trajectory (V c = V m ) or in the adiabatic case, eq. (pltjj) 
would lead to 

^ + m = Trfc, k = 1,2,3.... (47) 
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The violation of adiabaticity modifies this condition, and apparently, the bigger the differ- 
ence of the potentials V m and V c the stronger the deviation from (|47|). 

In the case of the Earth profile, the baselines L m and L c are correlated: they are 
determined by the zenith angle, Q u of neutrino trajectory: 

L m = R\ cos 0„| - L c /2, L c = 2^/R 2 cos 2 Q u -(R 2 -R 2 ). (48) 

Here R = 6370 km is the radius of the Earth and R c = 3486 km is the radius of the core. 
Using eqs. ([45)1. (f4~S|) we obtain that the condition (f4l)|) is satisfied for 

| cosQ u \ magic « 0.88, and | cosQ u \ magic « 0.98. (49) 

(For this estimation we took the potentials according to the average densities in the mantle 
and in the core p m = 4.5 g/cm 3 and p c = 11.5 g/cm 3 .) Thus, there are two magic trajectories 
for neutrinos that cross the core. They correspond to the baselines L = 2R\ cos W |, L magic = 
11210 km and L mag i C = 12485 km. For these baselines the total oscillations phases, c + 2</> m , 
equal ~ 4n and ~ 6tt. 

These new magic baseline could of interest for the atmospheric neutrino studies. 



6 More magic baselines 

The interference, and consequently, dependence on the CP-phase vanish also when Aa = 0. 
In this case the probability is determined by the solar amplitude: 

P(u e ^u e ) w \cos6 23 A s \ 2 . (50) 

Now the probability is determined by the solar parameters and can be used for measurements 
of 9 12 as well as and #23- 

In the constant density approximation, up to the phase factor, the atmospheric ampli- 
tude equals 

A a — sin 29^ sin (51) 

(The phase factor exp(— iAi 3 L/2) is omitted in (|5*Tj) . This factor should be restored if both 
solar and atmospheric amplitudes are non-zero.) So, the condition A a = gives the "magic" 
baseline (integer of the oscillation length in matter) 

Lmagic = <[(cOS20 13 " l u /kf + sin 2 29 13 ] ^ , (52) 

n — 1,2,... and here l v = 4irE/Am 2 3 . 

An interesting range of energies is below 3-4 GeV where the amplitude As is not sup- 
pressed too strongly by matter effects. Here, however, the "magic" baseline strongly depends 



9 



on neutrino energy. Therefore, a narrow energy neutrino beam should be employed or recon- 
struction of the neutrino energy should be done to suppress the interference and dependence 
of the probability on 5. For E u = 1 GeV the magic baselines are 1080 km (n = 1), 2070 km 
(n = 2), 3250 km (n = 3), etc.. For L = 3700 km that maximizes the solar amplitude, the 
"magic" energies are at 1.1 GeV, 1.5 GeV and 2.5 GeV. 



7 Conclusion 

The magic baseline (in the first approximation) is an integer of the refraction length. At 
high energies the latter approximately equals the oscillation length in matter. Therefore at 
the "magic" baseline the phase of oscillations driven by the solar mass splitting is 2tt, and 
consequently the solar amplitude vanishes in the transition probability. The interference 
of amplitudes, and consequently, dependence of probability on the CP-phase disappear. 
Defining the magic baseline as the distance that corresponds to vanishing solar amplitude, 
we have estimated various corrections to the equality L magic = l . The magic lengths in the 
non-uniform medium were discussed, and it is found that two additional magic trajectories 
exist for neutrinos crossing the core. We also discussed features of the magic baselines that 
correspond to zeros of the atmospheric amplitude. 
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